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Abstract. In games for which there exists a potential, the deviation-from-rationality dynamical model 
for which each agent's strategy adjustment follows the gradient of the potential along with a normally 
distributed random perturbation, is shown to equilibrate to a Gibbs measure. The standard Cournot model 
of an oligopoly is shown not to have a phase transition, as it is equivalent to a continuum version of the 
Curie- Weiss model. However, when there is increased local competition among agents, a phase transition will 
likely occur. If the oligopolistic competition has power-law falloff and there is increased local competition 
among agents, then the model has a rich phase diagram with an antiferromagnetic checkerboard state, striped 
states and maze-like states with varying widths, and finally a paramagnetic state. Such phases have economic 
implications as to how agents compete given various restrictions on how goods are distributed. The standard 
Cournot model corresponds to a uniform distribution of goods, whereas the power-law variations correspond 
to goods for which the distribution is more localized. 



1. Introduction 



Certain models in game theory [H 12 IHl UH 1201 12H1 1211 HOI HI] analyze the dynamics of decisions made 
by agents who adjust their decisions in the direction of higher payoffs, subject to random error and/or 
information (deductive "deviations-from-rationality" models, cf. [7j). These errors, which are essentially 
failures to choose the most optimal payoff, are understood to be intrinsic to the agents, and can be due 
to stochastic elements such as preference shocks, experimentation, or actual mistakes in judgment. In all 
of these contexts, the error is assumed to be due to intrinsic properties of the agents; i.e., the error is due 
to the agents making decisions that deviate from the true, optimal decision. These stochastic approaches 
commonly use a Langevin equation to deduce logit equilibrium measures for each individual agent. 

The logit measures had been earlier discovered by different reasoning [23 E]- Logit measures are 
exploited in [T1I2], where agents use a continuum of pure strategies with decisions perturbed by white noise, 
and additionally, agents can use a restricted form of mixed strategies which are absolutely continuous with 
respect to the Lebesgue measure. A purely dynamical approach is taken in j2Hl I2S1; where a continuum 
of pure strategies is assumed and the agents are restricted to using only these pure strategies with their 
decisions perturbed by white noise. If agents are restricted to play a single pure strategy at any given 
point in time as in [21112115 niixed strategies in mi2j can be viewed as a time-average of the agent's pure 
strategies. Since averages are effectively used in the stationarity equations of [H |2] , it is not surprising that 
the solutions for the logit measures are mean-field theoretic in nature insofar as they must find fixed points 
for these equations to solve them. A mean-field static equilibrium is also found for a discrete Ising-type 
logit model in 0. The dynamics of a random utility model are also analyzed in 0, and shown to have 
the same stationary state as with the mean-field analysis. The static equilibrium is that for a Curie- Weiss 
model. Such concepts will be examined below within the context of statistical mechanics. For the dynamics 
presented in this paper, we will assume that agents are myopic in pure strategy space; that is they play 
a single pure strategy at any point in time as in [23121]; and can only make infinitesimal adjustments in 
pure strategy space over infinitesimal time. The difference between the dynamical approach in this paper 
versus those in [H |21 HI HH 1201 123 121| is that a global approach is used \ in full exploitation of potential 
games [HSl- As a result, a single stationary measure results, instead of multiple coupled measures in [T1I2I. 
An explicit form of the measure is also obtained in contrast to the aforementioned references. 

It is noted that to allow agents to make effectively large moves in pure strategy over infinitesimal time 
(effectively what happens in P3I21) is akin to allowing mixed strategies, and as a result, Nash equilibria 
can be found much more quickly. As an example, consider pure-strategy quantities gi, q2, and which 
are positive real numbers. The dynamics of this paper does not allow a quick jump from qi to q2 since 
1^1 — 921 > 0, where | ■ | is the distance function on the line. However, with mixed strategies, an agent can 
shift from gi to (1 — e)gi + eq2 for small epsilon, and get a quick sense of the strategy q2 from qi and shift 
towards q2 in one step if it is better. This can happen since — [(1 — e)qi + eg2]|| < 2e, where || • || is the 
distance function on the Banach space of Radon measures. This is similar to classical, noiseless potential 
game theory when the Nash equilibrium is quickly and efficiently found to be a pure state (i.e., the global 
interior maximum of the potential function) since agents can make large unilateral changes in decision. 
This does not happen in the myopic, noisy dynamics presented here at nonzero temperature, as is evident 
from the resulting mixed strategy Gibbs/logit stationary state. The classical game-theory equilibrium (i.e., 
Nash equilibrium) is attained from the Gibbs measure at zero temperature. A thorough analysis of types 
of equilibria that occur in potential games can be found in [^ . 

In addition, an alternative view is proposed in this paper, in the context of potential games. The 
same Gibbs/logit equilibrium distribution that results from myopic, noisy dynamics can be derived using 



A single dynamical equation tracks all agents, instead of one equation for each agent. 
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certain axioms. This avoids explicit reference to error in agents' decisions and extends easily to games 
with discrete pure strategy spaces. These axioms are in fact those for equilibrium thermodynamics and are 
used to do large deviation theory. Therefore this alternative approach is empirical in nature: our model of 
agent behavior using thermodynamic axioms is based on many observations rather than the intrinsics (i.e. 
dynamics) of the game. It will be shown that this Gibbs approach to games with a continuum of strategies 
is equivalent to the aforementioned global, noisy, myopic deviations-from-rationality potential game. As in 
classical statistical mechanics, we will assume the observables (i.e., functions of strategies) to be C{Q), the 
continuous functions on pure strategy space. Hence the general mixed strategies will be the probability 
states in the dual of C{Q), Ri{Q) = {/i G C{Q)*\ /i > 0,/i(l) = 1}. That is to say the mixed strategies will 
be generalized to probability Radon measures on Q, which can include density functions with Lebesgue 
measure (such as the Gibbs measure), and delta functions for pure states (zero-temperature/pure-rational 
limits of Gibbs measures). 

Since we only use one dynamical equation, only one stationary logit measure is obtained and we will 
call this measure the Gibbs measure of the system, exactly as in statistical mechanics. This will also 
avoid confusion with the more general case presented in |2] of multiple logits for each agent. The 
Gibbs measure here is one for the entire system of agents in the given game which results from the global 
dynamical equation. This is a very interesting feature of potential games with noisy, myopic dynamics. 
It is stated in j^H] that "...in game theory, there is not a unique Hamiltonian, but rather each agent has 
her own Hamiltonian to minimize." This is true in general, but we'll show that we do have a unique 
Hamiltonian in the case of potential games. The idea of looking for something global to unify dynamics 
was touched upon in with the use of a "global cost function" in the 'minority game', but the technique 
is through the use of another dynamical method called "inductive reasoning" rather than errors to 
deductive reasoning. Within inductive reasoning, an explicit formula for the stationary measure in terms 
of the state of the agents does not appear. The equilibrium measures obtained here do have such an 
explicit formula: a Gibbs measure in terms of the potential. The potential was pointed out in to be 
the common part of the payoff functions that includes the benefits or costs that are determined by each 
agent's own decisions (e.g., the effort costs). Mathematically, the potential is just like the negative energy 
of all of the particles (i.e. agents) in a physical system (i.e. game) for a given state (i.e. point describing 
each agent's decision). 

The essence of the thermodynamic approach is that we are measuring many cases of the same game, 
as one would repeat an experiment many times to get average results; i.e., an ensemble average. The 
errors are not explicitly encoded at the (microscopic) level of individual agent decisions in this approach. 
However, this endogenous noise is implicitly assumed on a macroscopic level by requiring the ensemble 
average of the potential be constant. ^ The average values are due to making many measurements and 
averaging the results over the ensemble of the multitude of games being observed. Hence it's assumed 
that we don't know the details of probabilities of being in various states, but rather that we must deduce 
these probabilities from axioms of empirical observation. The size of the ensemble (i.e., the number of 
representative games we observe) must be very large so that statistical laws of large numbers hold, and 
that the empirical probabilities reflect the true intrinsic probability distribution of a single game. However, 
unlike the ideas proposed in [^12111211; we do not need to assume that there are a large number of agents 
and we do not treat a single agent as though she is a system in a 'heat bath' of agents, which implies 
an ensemble of agents and that we are only looking at one game. Rather, we treat the entire game as a 
system, and look at an ensemble of many identical games. This approach also allows an easy generalization 
to discrete/finite decisions. The partition function is introduced systematically and axiomatically, rather 
than as a post-hoc mechanism to solve a minimum free energy problem in inductive models as in I27j . 
The logit measure follows very naturally from the framework presented here. The limitation is that, for 



'If one game in the ensemble gains potential, another must lose potential which is counter to purely rational behavior. 
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now, a systematic introduction of a partition function (i.e., explicit stationary measure) is only justified 
for potential games, whereas ^1 121 123 is more general since agents have individual measures. 

An interesting aspect of this Gibbs approach is that the coefficient that determines the importance of 
errors in a stochastic potential game (the /i, = 2/(T^ in [H 121 or the Ijv^ in equation (jH))) is directly 
proportional to the usual notion of temperature in statistical mechanics. This was also discussed in |21] for 
individual agents. In fact, in the kinetic theory of gases in equilibrium, the temperature T is proportional 
to the average of the square of velocities v of the gas particles [T = k - (f ^) for a constant k). If a game has 
a potential that is the mathematical equivalent of such a gas, this equilibrium 'kinetic' theory states that 
the 'temperature' of that game is proportional to the average of the square of the temporal rates-of-change 
of decisions (T = c- {[dx / dt]'^) , where x G [x,x] represents an agent's decision). Hence larger instantaneous 
changes in decisions by agents correspond to higher temperatures of the game, and in general, temperature 
is a measure of the rate-of-change in decisions. In the dynamical sense Pfl |2] , temperature measures agents' 
deviation from rationality. 

A major application of statistical mechanics is to take infinite-agent limits from finite agent games, and 
an example of this will be done to show the classical result of oligopolies approaching perfect competition 
as the number of agents goes to infinity. This requires some modest technical restrictions to be made on the 
potential function so that the finite-agent Gibbs measures will converge in the infinite-agent limit (c.f. |42j): 
however these restrictions also make clear economic sense. A major advantage with the Gibbs approach for 
potential games is that, for various potentials, all of the rigorous techniques for phase transitions can be 
implemented as well as the process of renormalization of infinite-agent games. Renormalization would allow 
us to look at behavior of aggregate blocks of agents at various scales. Dynamical renormalization is also 
very recently being done in economics |T7||^, where large numbers of small firms are "agglomerated" and 
the dynamics of the resulting "meta-firms" are examined via a sandpile model. It will only be mentioned 
in this paper that renormalization can be done. The standard real-space renormalization is applicable to 
any potential game with the imposition of a lattice dimension, since all agents interact in many classic 
economic models. It would also be interesting to compare dynamical renormalization of the noisy myopic 
dynamics with the real-space renormalization of the Gibbs measure. Hopefully, the model shown here 
may be fruitful in understanding how local and global interactions between agents unify; that is to say a 
better understanding and more systematic treatment of economic interactions at small and large scales. 
Renormalization may link different economic scales as it does different scales in physics. 



2. Potential Games and Axioms 



Let us consider a game with a finite number of n agents, and all of these agents belong to the set A. 
At any moment in time, agent i E A can select an action or strategy Xi G Ai and the Xi are the decision 
variables. A configuration x of the system is any possible state of the system: x = (xi, X2, . . . , Xn), where 
each Xi G A^. The set of all possible configurations of the game is = ]^. Aj. A typical example of A^ is 
an interval of real numbers A^ = [x,x]. The development below is for continuum decisions, and discrete 
decisions will be accommodated later. 

In this article, we will only consider potential games [^. Let 

V{x) = V{xi, X2, . . . ,Xn) ■■ ^ ^ (-00, OO) (1) 

be the potential function for the game. Recall that a potential game is one such that every agent's payoff 
function Ui{x) satisfies 

^ = — (2) 

dxi dxi 
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for some function V{x). 



Certain conditions must be imposed upon V in order that probability distributions exist when the number 
of agents is taken to infinity in a hmit. These conditions are well known functional-analytic restrictions 
to ensure that infinite- volume limits exist (c.f. 02] )• Basic postulates for a game in equilibrium will be 
introduced (c.f. [SHI)- These are just the usual thermodynamics, and the probability of finding the game 
in a particular state can be derived from these postulates of empirical observation. 

Axiom 1 (conservation of energy). The game is isolated and agents act in such a way that the ensemble 
average of the potential is a constant. 

The 'isolated' requirement excludes outside sources from changing the potential for a given set of choices. 
Even though agents in a specific game in the ensemble try to maximize potential, the intrinsic noise will 
introduce errors into agents' decisions such that an increase in potential in one game will correspond to a 
decrease in the potential in one or more other games in the ensemble. Since purely rational agents attempt 
to maximize potential, this axiom precludes perfect rationality unless the ensemble average of the potential 
is the maximum potential (this is the zero-temperature case). 

Axiom 2 (equilibrium/stationarity). The isolated game is in equilibrium. This is to say that the probability 
of finding the game in any one state (i.e., a single representative from the ensemble) is independent of time. 

This is the usual commonsense notion of equilibrium. From Axiom |2l all macroscopic parameters (those 
that describe the system as a whole) are then also time- independent. Macroscopic parameters will be 
described in more detail later, but examples are the average output per agent, the average payoff, and the 
temperature. 

Axiom 3 (ergodicity). The isolated game in equilibrium is equally likely to be in any of its accessible 
states. 

An accessible state x is any configuration of decisions x & Q with V{x) = Vq, given a fixed observed 
potential Vq for the game (i.e., we look at all games in the ensemble with potentials of Vq). If a potential 
game with potential V has payoff functions Ui{x) for each agent i G A, then 



If the potential for n agents is restricted to a constant value Vq (i.e., we look at a constant-potential 
surface in fl), then the equilibrium dynamical system of P is ergodic. Fluctuations will cause a system 
to move among its accessible states over time with equal frequency, and ergodicity means essentially that 
the system, over time, will visit each point on the constant potential surface. As such, it is equivalent to 
assume all points on the constant Vq surface have the same probability of being observed. An interesting 
consequence of this that the condition dV/dt = does not necessarily imply an equilibrium (c.f. PP). If the 
system were at a local maximum of V (at say Xq) that isn't a global maximum, then Axiom El implies the 
system can't stay at the single point Xq. Rather, over time, it will move on the whole surface V = V{xo). 
Since V{xq) isn't a global maximum, the system will reach another point with the same potential but with 
nonzero gradient, and the dynamics can move the system towards a higher potential than V{xo). 

The approach to stochastic game theory in P produces a single Fokker-Planck equation that the joint 
probability distribution of decisions on Q satisfies. There, the Langevin equations for the joint distribution 
are 




(3) 




(4) 
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where the Wi are zero-mean, unit-variance normal random variables and z/ is a variance parameter. 



The following result very similar to [Ij. The major difference here is that we are looking at a dynamics for 
all agents who only use pure strategies, whereas in P single-agent dynamics are analyzed separately, and 
they use mixed strategies. The point of the following result is to show the stochastic dynamical approach 
in global form leads to the same Gibbs measure as that derived from Axiom ^ Axiom |2l and Axiom El 
The proof is in Appendix ^ Again, it is important to note that the dynamics below assume agents only 
use pure strategies and can only make small-distance moves over pure strategy space within a small time 
period (myopic agents). 

Proposition 4. Let f{x) be the joint density over decision space Q for a potential game with a finite 
number of agents n and potential V. Consider the dynamics 

dx = VVdt + iydw{t), (5) 

where x E Q, dx = {dxi, . . . , dxn), W = {dV/dxi, . . . , dV/dxn), and w{t) = {wi{t), . . . , Wn(t)) with the 
Wiit) standard Wiener (or white noise) processes which are identical and independent across agents and 
time. Furthermore, the Wi{t) have mean zero and variance one and reflecting boundary conditions ^ are 
used. Note that no conditional averages are done on V above, as opposed to [1112]. This indicates that only 
pure strategies are being played. 

If the process x{t) satisfies the dynamics of (01), then the joint density satisfies the Fokker-Planck equation 



df{x, t) j^^^ .^^^^^^ . ^ ^ VV(x, t) (6) 



and the corresponding equilibrium measure for variance is the Gibbs state 



In statistical mechanics, the term in the exponent of (jZl) is —E[x) / (kT), where k is Boltzmann's constant, 
T is temperature, and E{x) is the energy of configuration x. Hence the analogy of a potential game to 
statistical mechanics is that u"^ (deviation from rationality; influence of the noise in dynamics (@1)) is 
proportional to 'temperature' and the potential V is the negative 'energy' of the system. 

Now the Gibbs measure above can be derived using only the axioms of thermodynamics Axiom ^ 
Axiom 121 and Axiom El The proof is the standard proof of the 'canonical distribution' in any text on 
thermodynamics, and is the one that finds a state (i.e. probability measure) that maximizes entropy 
under the condition of a constant ensemble- average energy. The variational form of this maximization 
problem is the Gibbs variational principle of finding a state that maximizes the Helmholtz free energy. 
It is interesting to note that the problem of finding a state that maximizes the Liapunov function (for 
dynamical equilibrium) in P is identical to finding a state that minimizes the Helmholtz free energy of a 
given potential game. This is exactly what the Gibbs variational principle accomplishes. It is no surprise 
that the Liapunov function in [T] is the negative Helmholtz free energy for corresponding potential game 
with no explicit noise. This fact is stated here as a contrast to Proposition HI 

Proposition 5. Let f{x) be the joint density over decision space VL for a potential game with a finite 
number of agents n and a potential V. Suppose the three thermodynamic axioms AxiomU^ Axiom\^ 



'This requires zero time derivatives on the boundary, specifically that H51(l be satisfied for boundary points x. 
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and Axiom^^ hold. Then the equilibrium measure for this potential game is the Gibbs state at inverse- 
temperature (3 = (kT)^^: 

fix t) = fix) = "^^^^^^^^^ (8) 



Here the constant /5 in f{x) is determined so that the mean potential J^V{x)fix)dx = V, where V 
is the ensemble- average observed potential (i.e., the sum of the potentials of each game in the ensemble 
divided by the number of games g in the ensemble (Vi + ■ ■ ■ + Vg)/g). This is equivalent to the notion of 
considering one game to reside within the 'heat-bath' of the other games in the ensemble with a fixed total 
ensemble potential Vi + ■ ■ ■ + Vg. 

We consider equilibrium states of an infinite-agent game to be the measure (s) resulting from the appro- 
priate infinite- agent /volume limit of the Gibbs state. 

In the case of potentials that fit the usual notion of an 'interaction' ^ in statistical mechanics (c.f. ||42j for 
the appropriate machinery), the infinite-agent measures will be consistent with the finite-agent measures 
via the Dobrushin-Lanford-Ruelle equations. Such techniques can be used directly to analyze possible 
phase transitions for an infinite-agent potential games. We can then consider equilibrium states of an 
infinite-agent game to be the appropriate infinite-volume limit of the Gibbs measure. 



3. Example: Cournot Oligopoly, Local Demand, Perfect Competition and Collusion 



The simplest oligopoly model is due to Augustin Cournot (1838). There is one homogeneous good with 
demand function piQ), where Q is the total quantity of the good produced. Given an oligopoly of n firms, 
if each firm produces an amount qi of the good, then 

n 

1 

For the sake of later analysis, each firm's production q^ will be scaled between a minimum production g > 
and a maximum production q > q. We will assume each firm produces a sufficiently large number of the 
good so that the q^ can be regarded as continuum variables q < qi < q, I < i < n. For example, qi = I 
can represent the production of a sufficiently large number of goods. It is noted that smaller production 
would be handled in the discrete case as in Section HI 



Each firm uses quantity qi as its strategy based on the payoff function 

Il^ = q^piQ)-Qiq,), (10) 
where Ci is the ith firm's cost function. We will assume a linear (inverse) demand function 

p{Q) = a--Q, (11) 
n 

with constants a > and 6 > 0. Notice that b is divided by n so that demand is based on the average 
production. Thus demand stays nonnegative for large n. For example, if each firm were to produce q, 
p{nq) = a — bq is well-behaved and non-trivial (i.e., doesn't go to negative infinity or zero). 



It is noted that the potential for a Cournot ohgopoly is not such an 'interaction', and other techniques must be used. 
This is because of the 1/n factor; see Sectional for the exphcit potential. 
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Constant marginal costs (i.e., C- = c, for a constant c > 0) will also be assumed, so that Ci{qi) = cqi + di, 
and the payoff functions are now 



-b " 



= qi — ^ qj + (a - c)qi - di. (12) 



n 
1 



In the case of collusion, firms would try to maximize total industry profits 

n }) ^ n n 

n = ^n, = ^qiqj + {a- c)'^qi-'^di. (13) 

1 ^ i,j=i 1 1 

Both of these cases constitute potential games. A potential for the oligopoly in equation (fT^ is 

^o = -^Yl 1^^^ - 1^ E + - ^) E * (14) 

and a potential for collusion is the same as the collusion payoff function (except for the irrelevant constants 
di) in equation (fT^ . 

= ~^ E + - c) ^ (15) 

id 

These two potentials are two cases of the potential 

To get the potential in equation (UK), replace b and b in equation ((T]) with 6/2. To get the collusion 
potential in equation (fT^ . set 6 to zero. ^ 

The Nash equilibrium for the Cournot oligopoly with potential V above is 

a — c , , 

q* = . . 17 

' 2{b + b/n) ^ ' 

For non-collusion, q* = [n/ {n + 1)] [(a — c)/6] , and for collusion q* = (a — c)/ (26) are the Nash equilibria. ^ 
Note for collusion, the output is half of that for the non-collusive model, and hence firm profit (jl2p is less 
than it would be for the non-collusive equilibrium. These classical equilibrium solutions will be recovered 
in the perfectly rational case (zero temperature, /3 = 0) below, using the Gibbs approach. 



The partition function for n agents is 

h ^ h ^ . 

(18) 



^n = j W dqi exp 

i=i 



-13- QiQi - E + ^(^ - c) E * 



^It is interesting to note the effect of collusion pertaining to rationality. If a — c = in then collusion has the effect 
of doubling (3 in the partition function H18() . The collusive model is simply the non-collusive oligopolistic model at half the 
temperature 1/(2/3) = T/2. In the infinite-agent limit, the oligopolistic model is the same as perfect competition. Hence 
collusion is more "rational" behavior than perfect competition! 

^These equilibria look like those for total industry output Q = '^qi in the standard literature, but as described below 
, the scaling is at the level of firm production instead of industry production (where firm production would go to zero as 
n increases). 
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The infinite-agent free energy 



F{f3,a,b,b,c,q,q) = lim ln(Z„,/,). (19) 

- n-»oo (jn 



is calculated in Appendix ^ to be 

b{q + q)[b{q + q) - 2{a - c)] \n{q - q) 



F{f3,a,b, c,q,q) 



^ (20) 



— min —l3by'^ + ln ( sinhc 
p {/g(— 00,00) V 



m - q)y + (3^\^{b{q + q)-{a-c)} 



where sinhc(x) = sinh(x)/a:, sinhc(O) = 1, and it is noted that the free energy F is independent of 6, which 
can be set to zero. It is also shown in Appendix El that the expected value of any agent's output is 

, , d ^ a-c ^ ym{P,a,b,c,q,q) 

where < \{q + q)/2 — (a — c)/(26)|, ym < when (a — c)/{2b) > (g + '?)/2, and ym > ^ when 
(a — c)/(26) < (g + q)/2. Hence the Gibbs equilibrium is always between the average (g + g)/2 and the 
Nash Equilibrium (a — c)/(26). In the completely irrational limit 

limte) = (22) 

and the output is pushed to the average as would be expected for uniformly random behavior. In the 
completely rational limit, 

a - c 

.<^<«. (23) 



lim {qi) 



2b - - 26 

a — c _ 



which is the classical Nash equilibrium for an infinite number of agents in ()17|). 

From the properties of ym, we see that agents adjust output down (up) from g* = (a — c)/{2b) to- 
wards the average output qav = {q + q) /2 li q* itself is larger (smaller) than qav The magnitude of the 
adjustment increases as the deviation from rationality (i.e., temperature or 1//3) increases. This shows 
that a collusive-type behavior occurs insofar as agents produce less than the purely rational equilibrium 
output q* when q* > qav When q* < qav, agents will produce more output than g*, which is anti- 
collusive / supercompetitive. 

It is also shown in Appendix El that the volatility, which is x = P~^d'^F/dh'^, is zero at zero temperature 
(as expected for perfect rationahty) and goes to (g — g)^/12 as temperature goes to infinity (i.e., randomly 
acting agents have independent, uniformly distributed outputs qi). An agent's payoff can be computed 
using (g) and x- Since (11^) = lim„^oo g^l^V"^ Z]j=i (Ij + (« - c)] and (life) = (Ej), 1 < i < n, we have 

(Uk) = lim --X + {qk){h - b{qu)) = {qk){h - b{qk)). (24) 

n— >oo fi 

The maximum payoff occurs at the Nash equilibrium (g^) = g*. Furthermore, from the explicit formula for 
(gfc) it is evident that (11^) is an increasing function of /9. Payoff increases as the deviation from rationality 
(i.e. temperature) decreases. 

The model presented above is the standard Cournot model, which assumes the inverse demand is sym- 
metric among the agents. In reality, this is not typically the case. For example, water prices in parts of 
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a city with higher elevation may have a higher price because of additional costs of pumping the water. 
Likewise, competition among agents may not be perfectly uniform. An example of this would be firms who 
are more competitive against neighboring firms than with those far away, such as a local restaurant. It 
draws customers from mostly one area, and would be in greater competition with other nearby restaurants. 
The success or failure of distant restaurants would not have any affect. A single member of a restaurant 
franchise competes with local restaurants, but colludes with other more distant members of the parent 
company. In contrast to this would be a type of regional collusion, where all of the firms in any given 
locale may work in collusion with each other in an attempt to outcompete distant firms. 

To formalize these ideas, notice that an individual payoff function can be written as the collusion potential 
minus all other agents' payoff functions: 

nfc = K-5Zn,. (25) 

This shows an agent in oligopolistic competition will offset the collusive potential by subtracting the gains 
made by all other agents. A smaller degree of competition would correspond to subtracting fewer of the 
other agents payoffs. 

Example 1 (Local Collusion Within a Global Oligopoly). Suppose n agents are assumed to lie at integer 
points on the real line. If agents work collusively in disjoint groups of three, the payoffs are of the form 



--gi^i + (a - c)qi 



k = 3m, for integer m. (26) 



i=k,k±l i=k,k±l j=l 

A potential is Vic = Vo — (b/n) J2k=3mi=k±i where Vo is the oligopolistic potential in (fT^ . Although 
the potential is not translation invariant, if 'block outputs' Qk = {(lk~i + Ik + Q'fc+i)/3, k = 3m, are used 
along with periodic boundary conditions (g„ = qi) then a renormalized potential would be translation 
invariant. 

Example 2 (Local Oligopoly Within Global Collusion) . Suppose now that agents collude with non-neighboring 
agents and compete oligopolistically with neighboring agents within the same setting as Example ^ Here, 
however, the agents are competing oligopolistically with nearest neighbors, and not in disjoint groups of 
three. The payoffs are 

ni" = ^ n, = K - nfc_i - n^+i, (27) 

where we define IIq = and n„+i = to get correct payoffs. The potential is 

b 

Vio = Vc+-y2q^q^+i (28) 
1=1 

b^ , 6"-' 



n ^ — ' ^ — ' n 

1=1 



It is shown in Appendix Elthat the b term in the potential V of (|T6|l has no significance in the infinite-agent 
limit because the number of terms in the summand is of order n, and as such it can be set to zero. This 
is the case for the last nearest-neighbor sum in (j2EI), and in the infinite-agent limit, the free energy is the 
same as the collusion model. To maintain local oligopolistic competition in the infinite- agent limit, the 
interaction term b/n must be changed to a constant 6 > that does not vanish as n — > oo. 

Example 3 (True Local Oligopoly Within Global Collusion) . As mentioned in Example El an appropriate 
potential that maintains local competition in the infinite-agent limit is 

b n—l 

vjo = -- E + (« - c) E ?i + E (20) 
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where S > 0. 



The 5 term in ()30p only involves local interactions, that is, only terms within a fixed distance (of one) 
appear in the summand. Such terms in the potential are also called local interactions. 

Example 4 (Oligopoly With Stronger Local Competition). It is noted that in the infinite-agent limit, 
replacing b with 6/2 in ()30p is equivalent to changing the collusive potential component Vc to the standard 
oligopoly potential Vq. We can then interpret the (infinite- agent) free energy generated by the potential 
Vio to be oligopolistic competition with increased local competition. 



The potential in (I^Uj) is also generated by payoff functions 



n 



lo 



2b ^ 

n ^ 

j 



+ 6 



Qk-l + -7— Ik + Qk+l 

on 



(31) 



where go = Qn+i = applies. The dynamics generated by the in and the in ^Tl\ are identical, 
since they generate the same potential. Likewise, the infinite-agent free energies are the same, and as 
argued previously, the term bqk/{6n) will have no effect in the infinite-agent limit. 



What is interesting about 1)311) is that an inverse demand function appears: 

b 



Pkiqi 



PiQ) + 5 



Qk^ 



-1 + -^Qk + Qk+l 

on 



(32) 



where p{Q) is a standard linear inverse demand as in ()1H). The remaining term only involves terms local 
to agent k, and is called a local (inverse) demand function. Also notice that payoffs cannot be described 
by a single global demand p{Q), and each agent has her own demand function p^. 

The reason for the consideration of the addition of local demand as with ()30|) is that a phase transition is 
possible in contrast to the lack of one for the standard Cournot oligopoly. This will be shown in Section |3] 
for a discrete model which is a version of the Cournot model with local demand. 



Finally, we consider a situation where local demand for an agent's good is less affected by more distant 
firms. 

Example 5 (Power-Law Falloff of Local Demand) . Imagine agents residing at integer coordinates in a box 
in (i-dimensional space; call this set Q G Z,'^. Consider the local inverse demand functions 

Pfc(gi,...,gn) = -26 Yl |^5-|^ + («-c) + 2(5 (33) 

j'^Q-j¥=k ' i:\i-k\ = l 

where \i — j\ is the distance between agents i and j, and w > 1. Here, the (5-term again represents 
heightened local competition, but the 6-term is in a different form. The output of agents more distant 
from agent k has less of an effect on the inverse demand of agent k. This can be due to a variety of issues. 
Local agents may have more of a local market share because more distant agents could have to pay more 
for shipping, which can increase the price of their good in more distant regions. A variation of that is the 
good is primarily bought and not shipped, with customers preferring to travel less distance to purchase the 
good. This is in contrast to the previous term p{Q) in (|1H) which assumes goods are circulated uniformly 
in a global market, with no barriers that would put more emphasis on local output. 

A potential for payoffs 11^ = qkPk is 

= E + (« - c) 5^ g,, + 5 5^ g.g„ (34) 
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where indicates a sum over distinct pairs of sites in Q and indicates a sum over nearest neighbors 
= 1- 

It has been shown that a very interesting phase transition occurs for a discrete version of (jH^ above, 
which will be outlined in Section 0] This phase transition is very distinct from the phase transition 
speculated for the discrete version with the standard p{Q) term. We then see that an interesting feature 
of the Gibbs model of Cournot competition is that it distinguishes between models with uniform and 
nonuniform distribution of goods for certain ranges of parameters /3, b, and S. 



The inductive minority game (see, for example, OH] ) is a version of the El Farol bar problem 
01 127j that was introduced to analyze cooperative behavior in markets. It is introduced in the context of 
inductive reasoning and non-equilibrium statistical mechanics, but we will look at it in another way since 
the utility functions can be derived from a potential. An interesting point is that the dynamical equations 
in the inductive reasoning view do not satisfy detailed balance and there is a dynamical phase transition 
separating an ergodic from a non-ergodic process (c.f. the discussion in p!BJ). Similarly, there are Glauber 
dynamics for the Ising model which use sequential site updating and do not in general satisfy detailed 
balance, but nevertheless do converge to the appropriate Gibbs measure by stationarity (see |43| for an 
example). The lack of detailed balance usually suggest a study of dynamics, but the advantage of a Gibbs 
approach is that the dynamics yield an explicit stationary equilibrium state. 

In equilibrium analysis, many models exhibit a unique translation invariant infinite- volume state (hence 
ergodicity) at high temperature, while at lower temperatures there are multiple states and non-ergodic 
translation invariant limiting states exist (c.f. jl21 for details). As such it seems reasonable to consider 
a Gibbs model of the minority game. First, for the sake of background, we will introduce the inductive 
reasoning form of the game as presented in jTH]; and then contrast it with the Gibbsian version. 

The inductive game consists of n agents, labeled i = 1, . . . ,n. For simplicity, we can picture these agents 
lined up at integers on a one-dimensional axis, hence the agents are situated at sites A C Z. ^ At each 
discrete time step t, each agent i must make a decision cij G { — 1,1} (such as to 'buy' or 'sell'). The 
agents decide to buy or sell based on public information I{t) that is available to all agents at time t (c.f. 
the discussion in JH]). Note that at a given time, all agents receive the identical piece of information. 
Examples of such information could be the state of the market (decisions off all agents, the decisions of 
agents on odd sites, the block/regional- averages of decisions on blocks of length n/100, etc.), the weather 
forecasts, change in government regulations, etc. 

Profit, or utility, is made at time step t by agents according to the formula 



with 6 > 0. Because the coefficient —b < 0, the profit Ui{t) > when ai{t) = — sgn [^"^^ ^^(t)]. The 
justification is that in a market, if the majority is selling it is profitable to buy and vice-versa. Further 
notice the payoff function is scaled by 1/n. It is mentioned in jTH] to exist for mathematical convenience, 
but we will see later that it is needed for stability of the thermodynamic limit. Economically it makes 
sense to scale this way so that the inherent inverse demand function stays finite (see Section E)). 

Because of the form of the potential, this model is not a d-dimensional lattice model for any d. We have chosen to use 
d = 1 for visualization, and later, a true d-dimensional term will be added to the potential. 



4. Example: Gibbsian Minority Game 




(35) 



■t=i 
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Now the element of public information will be introduced. It is noted in (23] that a phase transition 
occurs without the element of information, and all that is needed is stochastic mixed strategies. In 
[T^ . information is randomly drawn at time t from a set {/i, . . . , Ip}. For example, at time t a number 
fi(t) G {1, . . . ,p} is drawn randomly and the information J^(f) is given to all agents, so that I{t) = I^{^t) 
with 1 < nit) < p. 

The information I{t) is converted into a configuration of n trading decisions o_{t) = {cri{t), . . . ,(T„(t)). 
These are determined using 'response strategies', which are look-up tables or functions IHj^s = {R]g, ■ ■ ■ , -Rf g) G 
{ — 1, 1}^. Each agent i has S decision making strategies ^Hj^g with s G {1, . . . , S}. The strategy that each 
agent i plays depends on the pubhc information I(t). If I(t) = I^{t), agent i will choose a strategy s and 

then lookup the decision -Rf^*^ = ±1 in the vector ^i^s- Agent i will then make the decision ai{t) = -Rfg*'', 
which depends on her choice of strategy and public information. 

Finally, the last ingredient is for agents to evaluate strategies based on the market performance of the 
strategies. At time t, strategy is profitable if and only if i?^^'' = — sgn [X]r=i '^ji^)]- Each agent i will 
then measure the cumulative performance of each strategy lHj,s, 1 < s < 5 by a profit-indicator pi^g which 
is updated at each time step t by 

p..(t+i)=p,.(t)-^i?i:f (36) 

j 

Note that at time t, all strategies are updated for each agent. Likewise note that the above is an evaluation 
scheme of strategies, and the strategies R are coupled to the actual decision ai(t) (since ai(t) appears in the 
sum). Additionally, these hypothetical decisions are being tracked in time, even though some aren't actually 
made. At time t, each agent i now chooses strategy Sj(t) which has the best cumulative performance Pi,s(t) 
up to that time. 

Given the external information n{t), the dynamics of the minority game are given by 

p,,,(t + l)=K,(t)-^i?g*)/l(t), 

3 

Si{t) = argmax^g|i_ .5.|j9i_5(t). 



Some interesting observations about the above approach need be made. Firstly, a phase transition 
occurs in the above inductive model even though the model is a totally frustrated one. The fact that 
something stochastic is what introduces the phase transition makes sense since this is what happens in the 
setting of spin glasses and replica symmetry breaking (c.f. [121 123 E] ) • It has also been shown in |221EZ1 
that the inductive reasoning model seems insensitive to the types of long range interactions that may 
occur. Speculatively, it seems that the global nature of public information can outweigh the interactions 
among agents, just as a global magnetic field can alter the susceptibility/volatility in a fully frustrated 
model. However, there are some fundamental physical and economical differences between the long-range 
interactions presented in [221 The economical differences were described in Section |2] as local versus 
global inverse demand and/or availability of goods. The physical differences between such interactions are 
significant, as will be shown below. 

This motivates us to study the minority game in a different light than the inductive approach above. Now 
we will introduce the Gibbs model of the minority based on statistical observation of deductive reasoning 
and show how this model distinguishes local and global demand. The Gibbs approach is not inductive: 
only the values of the agents' decisions are considered. Which inductive strategy an agent may be using is 
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irrelevant. In other words, we do not distinguish inductive strategies 9^*'® that an agent may be considering, 
but only look at her final decision (Xj to buy or sell. 

In the simplest Gibbs model, there is no phase transition, but the infinite-volume volatility (c.f. Ap- 
pendix ^ for the continuum version) does display nonrandom behavior at finite temperature. There are 
phase transitions in the Gibbs model if local interactions among agents are added which represent greater 
competition among local firms. This will be pursued below. 

Notice that the utihty functions in can be deduced in the sense of ^ from the potential 



where 6 > and there are n agents. The agents can be pictured to lie on integer points of some lattice, 
but dimensionality of the lattice is irrelevant since all agents interact. If short-range perturbations are 
added, the dimensionality of the lattice is then relevant. In the basic model presented here, the volatility 
is differentiable in the global field (i.e., no phase transition; see below). However, we will momentarily 
digress from the basic model in (j38|) and look at variations in an attempt to determine what may cause a 
phase transition. 

Consider this new model (in (jHHj) below) to lie in two dimensional space, and the addition of nearest 
neighbor ferromagnetic /aligning interactions. These local aligning interactions correspond to the (local) 
inverse demand for an agent 



where 6 > and 6 > 0. When n is large, the local aligning interactions having coefficient 6 correspond to 
the ith agent's inverse demand function increasing in response to local supply increase (i.e. selling, or the aj 
changing from +1 to —1 for |j — «| = 1). On the surface, this seems contrary to usual supply and demand, 
since inverse demand should be a decreasing function of supply variables aj. ® Mathematically, these 
interactions give more probability weight in the Gibbs measure to configurations with more alignments 
between local agents. This could be interpreted as a reward when nearest agents to do the same thing 
(i.e., align: all buy or all sell). It could be a purely psychological reward for "following the herd", and 
not part of the monetary economic demand function, since utility or payoff functions need not be based 
on monetary payoff alone. It could also be a form of monetary influence where agents are willing to pay 
other agents to go along with their ideas. Since an agent only has a relatively small amount of money, they 
would only be able to payoff several agents. However, it is important to note that the flrst term of (jHI)]) 
may penalize the agent for alignment. Ultimately, whether or not an agent beneflts from aligning depends 
on the phase of the system. 

A potential for the model with payoffs in ()39|) is 



where b > 0, 6 > 0, i, j E Z^, and the last sum is over nearest neighbors \i — j\ = 1. It has been speculated 
that the model in !^40\i has a transition as T = 1/(3 ^ from a fully frustrated state to a ferromagnetic, 
or fully aligned, state provided b < AS ^U]. If 6 > 45 there would likely be no transition and only full 
frustration as with (jHHj) [TU]. Completely aligned states should occur in some cases, and there should not 
always be only fully frustrated states. 

^One interpretation in Section 13 is that such ahgning terms are due to heightened local competition. Aligning terms result 
from agents subtracting the utility/payofF of nearest-neighbor agents from their own utility functions. 




(38) 




(39) 




(40) 
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In the case of overhyped introductory price offerings, essentially everyone is in a frenzy to buy and we 
would see the configuration of cXi = —1 for all i. If news got out that a company lost all of its worth, 
everyone would sell and we would see the configuration of (Xj = +1 for all i. The aligning interactions are 
the only way to support such states without the addition of a field term {a — c)^ai, as the model in 
has no phase transition (cf. Appendix IH]) . The addition of such a 'field' term, which has demand and cost 
components, will support such states in the absence of a phase transition. The demand/cost component 
should be able to support such states, since such extreme market reactions result in pricing adjustments 
(changes in a). However, it was shown in Section El that the aligning terms can be interpreted as increased 
local competition. Such effects can become large in such a frenzied environment, and large alignments 
may occur no matter how small the demand/cost term a — c may be if the temperature is low enough (i.e., 
agents are not deviating much from "rational" behavior). 

The model in (j4(Jj) is similar to the models studied and summarized in pUllTTllHHj . The salient difference 
is that innHHlESI study a long-range 'Coulomb' {w = 1) or 'dipole' {w = 3) on the two-dimensional square 
lattice 



where the i,j G are two-component vectors with integer coordinates. These models are very rich, and 
interesting in an economic sense: the local demand functions should depend more strongly on local agents, 
and not depend significantly on distant agents. This suggests that output of goods local to an agent has 
more of an affect on that agent than goods that are distant from the agent, i.e., that local inverse demand 
(cf. Section depends more on local supply. The model in PUI) assumes that local inverse demand is 
identical for each agent, and that distant goods are equally accessible to an agent as local goods. There is 
a fundamental difference between these interactions. An important difference in two dimensions between 
pn|) above and models with falloff potentials where n is replaced with \i — jI"", w = 1,3, is that the latter 
models can have a phase transition to an antiferromagnetic phase at low temperatures and the former 
model has no such phase transition and only totally frustrated ground states. Antiferromagnetic phases 
correspond to a checkerboard of +l's and — I's in the plane, which is ordered behavior of agents. The 
agents with potential (j4Up have no such ordered behavior: +l's and —I's are seen equally but in a random 
unordered fashion, which reflects that supply is equally accessible to everyone. 

There are four types of phases in the dipolar model (cf. for details and references). The first two 
occur below a critical temperature Tc(A) that depends on A = 6/b. When A is smaller than A^ ^ 0.425, 
the ground state is the antiferromagnetic one, consisting of a checkerboard of +1 and —1 values of the a. 
When A > A^, the transition is from the antiferromagnetic state to a state consisting of antialigned stripes 
of width one, representing large areas of buyers and large areas of sellers that persist in some direction. 
As A gets larger, the stripes grow in width. For example, stripes of widths two, three, and four appear at 
Afc ^ 1.26, 2.2, 2.8, respectively. 

The other two phases occur when T > Tc(A). There is a transition to a disordered tetragonal phase 
above and near the Tc(A), and this phase is not paramagnetic. There are extended ferromagnetic domains 
characterized by predominantly square corners, and a typical configuration looks like a maze of +1 trails 
and —1 walls. For larger temperatures, the transition is to a paramagnetic phase. Similar phase behavior 
occurs for the w = 1 Coloumb model in ()41|) [33] • 

The antiferromagnetic phase is one with buyers and sellers in a checkerboard pattern. This hap- 
pens when the demand p{Q) is mostly due to nearest-neighbor oligopolistic competition; i.e., the terms 
~^^r\i-k\=i1k1i dominate agent fc's demand function. When there are stripes there are large areas per- 
sisting in some direction that are buying and other such areas that are selling. In this case, agent k has 
localized oligopolistic competition in one direction and heightened local competition in the other direction. 




(41) 
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As the stripes widen, there is heightened local competition within the stripes, and a localized oligopolistic 
competition among the groups of agents in the stripes. In the tetragonal phase, there are large unordered 
areas of buyers and sellers that look like a maze. In the paramagnetic phase, buyers and sellers tend 
to appear at random, and this corresponds to the standard Cournot model of nonlocalized oligopolistic 
competition. 

It seems intuitive that (jiUI) would not display stripes, but rather be fully frustrated or completely aligned; 
that is, if supply is globally accessible and local prices are essentially the same for all agents (such as a 
homogeneous market with no trade barriers and uniform distribution of goods) then there should not be 
large clumps of buyers and sellers. On the other hand, if there are no local interactions as in (jHHj) . then 
there should only be antiferromagnetic states or fully frustrated states with high probability, and no large 
clumps of buyers or sellers. 

The numerical evidence in l?7j for the inductive minority game (j37j) suggests it is insensitive to the 
types of long range interactions; specifically that the types in (j38|) of equal interaction strength show the 
same behavior as the fallofF types in (j4H) with 6 = 0. This is in contrast to the deductive model presented 
here. 

The stochastic nature of public information in the inductive model does induce a phase transition. It 
would be interesting to determine if the deductive model presented here with a random field component, 
which could be represented by random marginal cost fluctuations (cf. Section El), will yield a phase 
transition. Such random costs could, for example, represent the fees, which vary from agent to agent, 
charged to each agent per buy/sell. In jSHI; the addition of demand into a model of stock trade is analyzed 
by introducing an effective field as mentioned above. 

Evidence is found in j^Z] that the dynamics of the inductive minority game is intrinsic to resource 
allocation. The view of the minority game as a potential game makes clear why this is the case, in the 
sense that all potential games are isomorphic to congestion problems [^3] , and the potential for the minority 
game is a specific, discrete case of a supply-and-demand potential (cf. Section E)). 



5. Conclusion 

We have seen that potential games (with nice potentials), when coupled with deviations-from-rationality, 
result in models for which equilibrium values of all relevant (macroscopic) quantities such as output per 
agent, payoff per agent, etc., are determined by averaging over a Gibbs measure with the given potential at 
inverse temperature f3. The introduction of temperature in economics simply corresponds to the assumption 
that players deviate from the purely rational behavior of producing output in the direction of maximum 
increasing payoff. 

A very nice feature of the Gibbs approach is that it is accessible and demonstrates different patterns 
of behavior among agents for different types of distribution of goods. This is desirable when modeling 
trade barriers, shipping costs, and real-life models such as restaurants, where demand is localized. It can 
also be used in conjunction with renormalization techniques, which could model large-scale behavior of 
'agglomerated' agents, analogous to block spins for physical systems. 

It would also be of interest to pursue more depth insofar as the economic implications of different types 
of phases for the different models presented here. Since the author is not an economist, such an endeavor 
will be left to the experts. 
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Appendix A. Proof of Vector Fokker-Planck Equation 



Here we will show the derivation of a vector Fokker-Planck equation for the joint distribution given the 
Langevin dynamics in p. This proof is just a slight modification of the proof in We only consider a 
finite number of agents ?7, < cxd. 

In a potential game, the Ito equation for dynamics is dx = Wdt + adw{t) as in (0)). For a small change 
in time At, the dynamics can be written 

Ax{t) = x{t + At) - x{t) = VVAt + aAwit) + o(At), (42) 

with aAw{t) a normal random variable of mean zero and variance a'^At, and d{At) is the usual vector 
'little order' notation. ^ The configuration of decisions x{t) at time t is hence a random variable with 
time-dependent joint density f{x, t). Let h{x) : Q (—00, 00) be an arbitrary twice differentiable function 
such that h and Vh vanish on the boundary dQ of decision space Q. Note h does not explicitly depend on 
time. Then the expected value of h{x) (over phase space) at time t + At is 

E [h{x{t + At))] = / h{x)f{x, t + At)dx, (43) 
Jn 

where dx is the product measure dxidx2 ■ ■ ■ dxn on VL. Note above that x{t + At) on the left side represents 
the random variable over which the expected value E is to be taken whereas the x's on the right side 
are the points in VL. Equation (j42|) can be used to obtain another expression for (j43|) 

E [h{x{t + At))] = E [h{x{t) + Ax(t))] = E \h{x{t) + VV{x{t))At + aAwit))] + o(At). (44) 



The expression on the right hand side of (j45|) will be subtracted from the Taylor expansion of the right 
side of (HH) to yield the Fokker-Planck equation. We will proceed as follows. 



Let g{y) be the (joint) density of aAw{t), which is an n-dimensional normal density with mean zero and 
variance cx^At. The left side of (j44p is in terms of the random variables x{t + At), and the right side is in 
terms of the variables x{t) and y with respective densities f{x,t) and g{y)- Accordingly, the right side of 
(j44|) can be written as an integral over these densities: 

E [h{x{t + At))] = / h (x{t) + VV{x{t))At + ay) f{x, t)g{y)dxdy, (45) 

where M = (—00, 00). 

The Taylor expansion of the right side of (j45p yields 

h{x{t)) + V/i(f (t)) ■ [VV(f(t))At + 

1 . . . ^''^ 

+ - [W{x{t))At + y\ ■ D^h{x{t)) ■ [VV{x{t))At + y\'^ + d{At) \f{x,t)g{w)dxdy, 



R" JQ, 



where D'^h is the second-derivative matrix for h (i.e., the ij entry is d'^h/dxidxj) and the transpose of a 
row vector v is the column vector v'^. 



^v{At) is o{At) if liniAt^o Vi{A.t) / (At) — for each component Vi of v. 

^'^Rigorously, the map x{t) — > x{t + At) generates a differentiable global (stochastic) flow on decision space for nice 
potentials, and the global flow has a generator which will be shown below. 

^^For simplicity, we use the formal notation of differential stochastic calculus. 
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Integrating ()46p over the y eliminates terms linear in yi as well as the mixed second-order terms since 
the yi have mean zero and are independent. The expected value of wf is cr^At, hence fHUj) reduces to 



h{x{t))f{x,t)dx + At / \/h{x{t)) ■WV{x{t))f{x,t)dx 
n Jn ^^^^ 



2 f 

+ At— / V^h{x{t))f{x,t)dx + o{At), 
2 Jq 



where o{At) is the one-component version of o(At) and is the Laplacian. Thus the generator for the 
Langevin dynamics x{t + to) = ^'^x{to) is L = (VV") ■ V + (cr^/2)V^, where L acts on twice-differentiable 
functions on Q that vanish on the boundary. It is easy to show L is symmetric and the domain of its 
adjoint is the Sobolev space H'^{Q) for a well-behaved potential V. We view its adjoint L as a densely 
defined unbounded operator in L'^{Q,dx). The idea is apparent from the following. 

Using Green's identities, the integrals containing the terms V/i and V^h can be integrated by parts, and 
the boundary integrals vanish on the boundary dQ (since h and V/i do) leaving a new expression for ()47|1 : 



h{x{t))f{x,t)dx- At / h{x{t))W ■ [VV{x{t))f{x,t)]dx 
n Jn 



2 

+ Aty / h{x{t))V^f{x,t)dx + o{At). 



n 



This shows that the adjoint of L operates as L*(/) = -V ■ (/VV) + {a'^/2)V^f. Equations ^ and (gHI) 
are identical, hence subtracting them yields 

/ hixmfi^,t + At)-fix,t)]dx = At [ h{m)[-^-NVixit))fix,t)] + ^V'f{x,t)]dx + o{At). 
Jn Jn ^ 

(49) 

Dividing ()49p by At, taking the limit At ^ results in the Fokker-Planck equation (jH)) since the set of all 
such /I's described are a separating set for probability densities on (i.e., these /I's can weakly approach 
delta functions on fi). 

Let h[x) be an observable (i.e. function) on that is differentiable, vanishes on Sfi, and doesn't depend 
explicitly on time. For a stationary state /, the average of h[x) is constant in time: 

= h{x)f{x,t)dx 

= yMx)V- |-[W(f(t))/(f,t)] + yV/(x,t)| (50) 

= - jvh{x) ■ |-[W(f(t))/(x,t)] + y V/(f,t)| . 
Since h was arbitrary and V/i can be made to approach delta functions, for a stationary state 



-[Vr(f (t))/(f , t)] + y V/(f , t) = 0. (51) 

We then see that the Gibbs state in ((Tj) is the unique solution to (HH), and is the equilibrium state for the 
potential game. The negative of the Helmholtz free energy is in fact a Liapunov function for the dynamics 
and it can be used to show (see jT]) that the unique finite-agent solution equation ((Tj) of equation (j3Tj) is 
attained in the long run. 



19 



Appendix B. Solution of the Cournot Free Energy 



The free energy for the Cournot ohgopoly and collusion is obtained from the partition function in 
equation (fTH|l 



ZA b 



The limit of the free energy 



exp 



b b 



j,k=i 



will be shown below. First note that the limit does not depend on the value of b, since 

bq'^ /-\ bq^ 



n 



n 



(52) 



(53) 



(54) 



where = -F'„(0). It then suffices to set b to zero in equation (j^ . which will be done below to simplify 
a limit. For further manipulation, changing variables to g = g — 7, 7 = (g + g)/2, /i = a — c, results in 



Zn = exp 



n 



b h 



-pbYn - PbY + Ph-fn + \ jb + 



X 



exp 



-Q/2,Q/2]- 



n 2 



b h 

67 + -7 - o 
n 2 



(55) 



\{dqj, 



where Q = q — q. 



The solution can be found using a saddle point method, and the key is to use the identity 

«p2 _ 1 



ny/4p 



exp ( + ^^'^ ) '^^ 



(56) 



as with the simplest case in j^I]. Using the above in and setting 77 = 67 + 67/72 — /i/2 results in 



Zn 



exp 



46 \ bn 6 



exp 













{ I exp 




exp 






yj-Q/2 


n 









(57) 



where the order of integration was changed via Fubini's theorem. Consider part of the integrand above 



9n\ 



1 

QJ- 



Q/2 



exp 



Q/2 



n 



exp 



iq\l-t 

n 



dq = sine 



Qlb 



2 V n 



(58) 



where sinc(x) = sin(x)/x (sinc(O) = 1), 6 was explicitly set to zero, and a normalizing factor was added so 
that gtiO) = 1- Changing variables to t = tj ^Jn^ the partition function can then be written 



Zn = exp [-p'y{ky - h)n] 



X 



exp 



1 



f 2P 



sinc^^V^t^ 



(59) 



dt. 



Now a basic fact is needed to determine the limit of the free energy. 
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Lemma 6. Let f{x) be a bounded, continuous, real-valued, integrable function on the real line M and 
let G be an open set. Suppose there is a point xq G G for which \f{x)\ < f{xo) for all x ^ xq and 
limsup |/(x)| < /(xo) for all x outside of some neighborhood of Xo contained in G. Then 

lim / [fix)r^' dx = /(xo) (60) 



m—^oo 



G 



for any fixed integer k. Furthermore, for a function g{x) that is continuous at xq with gf^ G , 



m— >oo 

-oo 



Proof. First, by factoring ||/||oo out of the integral in ()fi()|l . we may assume /(xq) = 1 by considering the 
function //||/||oo- Choose a number cq > so that hmsup |/(x)| < l—2eo for all x outside of a neighborhood 
AT c G of Xo and / > on Af. Let = {x G G|/(x) > 1 - eo}, and J„ = (/(.[/ (2;)]'"+'^ rfx)i/". Since 

Im < f{xoY"''^'^~^^^"^{Ji^f'^{x) dxY^"^ for large m, limsup^^o^/m < /(xo). The inequality 



{i^r > (1 - e)' 



-m+fe-2 



/ f\x)dx-{l-2e,r+^-^ ! , 

J5i_E JG-M 



/^(x) dx 



(62) 



holds for all e < eo- Upon factoring out (1 — 2^ jg ^^^^ ^^iSit liminfm^oo > fi^o)- To prove (|UT|). 

the following is needed: 



lim 

m— >oo 



Ir-s,.9{x) [/(x)]"" dx 



lim 



L-s,.Jf(^T''dx 

°° \iG[f(^)r dx\ 



0, 



(63) 



where the last equality follows from the root test for convergence since the hypothesis implies the limit of 
the ratio on the right side of (jHSj) to the power 1/m is less than one. With this, reduces to 

js,_jix)[f{x)rdx 
iu[fi^)rdx ■ ^ ^ 

By the continuity of g and (jU^ again with G = Si^^, it is seen that holds. □ 



To apply Lemma IHl to the partition function in ()59p . the integration path along the real line must 
be deformed in the complex plane to a path on which the integrand is real-valued. As such, consider 
the function f{z) = exp [— /c(z + i?7o)^] sinc(z), where rjo = —PQrj and 1/k = pbQ"^ (changing variables 
t = QVbi/2 in (jKHjl : note r/o is increasing in h). 

Let z = x + iy, and f{z) = 7l{x, y) + iT{x, y), where x, y, TZ, X are real- valued. Then we are interested in 
the paths (x(s), y{s)) in the complex plane C on which X = 0. Since / is holomorphic and not identically 
zero, the gradient of X, VX, can only vanish on a closed set of isolated points in C. Thus X~^(0) consists 
of various differentiable paths in C. The imaginary part can be expanded 

X(x,y) = — { sin(—2/cx(?/ + r7o))k sin X cosh ?/ + ?/ cos xsinht/l 

+ cos{—2kx{y + i]o))[—y sin x cosh ?/ + x cos x sinh y] } . 

Lemma 7. There is a piecewise differentiable path z{s) such that I{z{s)) = 0. In addition, z{s) can be 
written {x,yz{x)) for small and large \x\, with lim^-^-i-oo ^2(2;) = Ci and lim^-^o ^2(2;) = C2 for constants 
Gi and G2. 
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Proof. Consider fundamental domains in the 2;-plane of w = sinc(z). The fundamental domains Dj, Dj 
shown in figure ^ are then separated by curves satisfying y = and tanx/x = tanhy/y, with x ^ 0. If 
X > 0, the implicit function theorem implies the latter curves can be written x = Cj{y) with Cj different iable 
and jvr < Cj{y) < (2j + l)7r/2, j = 1, 2, . . . . Let Co{y) = and for x < 0, the curves are reflected: c_j = —Cj. 
The branch cuts in the w-planes are then the images of C2j+i. The logarithm ln(sinc2;) will now map two 
interior points in the z plane to a complex number since sinc(— z) = sinc(2;). Branch cuts at |x| > are 
added in the z-plane on the x-axis to avoid discontinuities in a. The fundamental domains Dj lie below the 
X-axis between the curves C2j-i and C2j+i, whereas the D'j lie above the x-axis between c_2j-i and c_2j+i. 
The shaded regions of the Dj and Dj in the figure below correspond to the lower half of the tu-planes. 

Let a = a{x,y) be the angle of sinc(2) (i.e., the imaginary part of ln(sinc2;)). Figure ^ shows that in 
the lower half of the z-plane, a increases in x and is of order x. 




Figure 1. In(sinc2;) values 



Differentiating the identity for tan < 

da 2y (sin^ x + sinh^ y) — (x^ + y"^) sinh 2y 
dx 2 (x^ + 1/2) (sin^ x + sinh^ y) 



(66) 



da — 2x (sin^ x + sinh^ + (x^ + y"^) sin 2x ^^^^ 
dy 2 (x^ + y"^) (sin^ x + sinh^ y) 



and it is seen that If > (< 0) when y < (> 0). 
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The part in braces in ()65|) can be written 



a/x^ + Y sin^ X + sinh^ y sm{—2kx{y + r^o) + (68) 

and this is zero if and only if X is zero. For the case 7]q = 0, the integral in (j59|) can be evaluated using 
Lemma ini with the maximum of the integrand occuring at zero. In the case r]o 7^ 0, the paths in C on 
which f{z) is real- valued is precisely where 

—2kx{y + riQ) + a = m'7i, (69) 

where m is an arbitrary integer. 

Since y + fjo = a/{2kx) — rm[/{2kx) and a is of order x, it is evident that |?/| — *• 00 as x — *• on such 
paths unless m = 0. This is not surprising, since large \z\ behavior of / is dominated by the term in the 
exponential (thus giving hyperbolae-type curves). We will single out the m = curve z{s) = {x{s),y{s)) 
satisfying 

y + rio = a{x,y)/{2kx) (70) 
for reevaluating the integral in (fS^. Since a is of order x, the ?/-component of z{s), y{s), is bounded 
and has finite limits as x ^ ±00. The location of the curve depends on rjo- To see this, note that 
a{x,y) < Dx for a constant D with \D\ < 1, and the sign of D is negative in the upper-half and positive 
in the lower-half of the 2;-plane. If r^o > and the curve z{s) were in the upper-half of the 2;-plane, then 
y = a/{2kx) — tiq < D/{2k) — 779 < 0, which is a contradiction. Hence if r^o > (< 0), then the curve z{s) 
is in the lower (upper) half of the z-plane. 

The ?/-intercept of z{s) can be found using lim3;.^o+ « = and limaj^o sin(a)/Q; = 1. Near x = 0, z{s) 
can be written z = {x,yz{x)), from which 



a sin X cosh + X cos X sinh y 

yJyOj + rjo = lim — — = lim 



x^o2kx x^o 2A;xa/x^ + UzV sin^ x cosh^ yz + cos^ x sinh^ y^ 
y^O) coshyj;0) + sinhy^O) 



(71) 



2%^0)sinhy^0) 

The solutions of the above equation in yz{0) are the same as the solutions in y to the equation y/[l — 
2k{y + r]o)y] = tanh(?/), except for an extraneous root of y = in the latter (since r/o 7^ 0). The solution 
for the y-intercept is unique, and this will be shown later. □ 



With the above, the integral in can be deformed onto the curve z{s). To apply Lemma IHl it is 
necessary to find the maximum of |/| on the full curve where 1 = 0, and to also show it occurs at a single 
point Zm such that f{zm) > 0. 

Equivalently, it will be easier to find the maximum of 

l/P = exp [— 2A;x^ + 2k{y + r^o)^] (si^^ x cosh^ y + cos^ x sinh^ y) / (x^ + y'^) (72) 

on the full X = curve z{s) = {x{s),y{s)). Since lim^^-i-oo f{z{s)) = lima-^-too f{x, yz{x)) = 0, a maximum 
occurs and is greater than zero. On the curve near this maximum, say at s^, the argument of / is constant 
(either Arg(/) = or Arg(/) = n) hence tan(Arg(/)) = X/TZ is constant. Since tan(Arg[/(x(s), ?/(s))])' = 
0, the tangent vector v{s) to the curve {x{s),y{s)) is orthogonal to V{X/7l) near Sm- 

Any maximum point of |/p on the curve z{s) must satisfy 

„_ d ,_ exp[-2x^ + 2(|/ + r7o)^] 

ax'^' (x2 + y2)2 ^73) 

X {— 2x(2A;x^ -|- 2ky'^ + 1) (sin^ x cosh^ j/ -|- cos^ x sinh^ + (x^ + y"^) sin2x} 
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and 



= A|/|2 _exp[-2x2 + 2{y + r/o)^ 



dy''' (x2 + 1/2)2 ^74) 

X {^{Ak{y + r]o){x^ + y^) — 2y) (sin^ x cosh^ ?/ + cos^ x sinh^ + (a;^ + y^) sinh 2?/} . 

Lemma 8. The conditions in (ffH|) and fl74p imply that the maximum of |/p on i/ie curve z{s) occurs at 
(0,^/2(0)), i/ie y -intercept of the curve z{s). 



Proof. Without loss of generality it can be assumed that ^ 0, since the integral in ()59|) is an even 
function of rj^. 

For the case rjQ = 0, the curve z{s) is the a;-axis and \ f\'^{z{s)) = e~^*^^^ sinc^x. Clearly (0, 0) is a global 
maximum of |/p on z{s). Hereon, the case r/o > is considered. 

We will consider the case x 0, and x = will be dealt with subsequently. Then (f7H|) is equivalent to 

2kx'^ + 2ky'^ + 1 . o 1 2 x x 

1 > sine 2a; = (sm x + smh y) > 0. (75) 

x^ + y^ 



The inequality sine 2a; > implies x G [— 7r/2, 7r/2] or x G ±[m7r, (m + 1)71 /2] for positive integers m. In 
order for (j75|) to hold, sin^ x + sinh^ y < sinc(2x)(x^ + y^), which is impossible for |x| < 7r/2. Any extreme 
points of l/p for which x 7^ must therefore satisfy |x| > vr. 

From the above argument, maximum points of |/p in the set So = {{x,y) : |x| < 7r/2} must satisfy 
X = 0. Any maximum point (x, y) ^ So will then satisfy |x| > tt. It is seen from fl70|) that ?/ + r^o < when 
X > and y < 0, which along with (ff^ . implies that \y\ > (tt^ + y^) tanh \y\ at an extreme point. That 
this is impossible shows the largest value of |/p on the curve z{s) will occur at x = 0. 

It only remains to show the intersection of z{s) with the y-axis is unique and is also a maximum 
point of I /p. Note that (|7H|l is satisfied when x = 0. Also note that (f7I|l reduces to the condition 
[4:k{y + rio)y'^ — 2y] sinh^ y + y"^ sinh 2?/ = As described by (fTTj) . ?/ = is an extraneous root since r/o > 0, 
and (fflj) is satisfied if and only if 

^ = ^ (- COthy + }^J+Q{h^-^^ giy). (76) 

As was mentioned in Lemma any solution y^ must be negative since t]o > 0. The function g{y) on 
the right side of (fTBj) decreases from bQ^/2 + (^(67 — /i/2) to Q{b'y — h/2) < as y increases from —00 
to zero. Therefore g{y) intersects h{y) = y/j3 a.t exactly one point yo < 0. This shows that there is a 
unique solution |/o < to (ffH|) . Consequently, the global maximum of |/p on z{s) occurs at {0,yz{0)), the 
y-intercept of the curve z{s). □ 



Combining the result of Lemma |H1 with Lemma El it is seen that 

F{/3, a, 6, c, g, q) = lim ln(Z„,^,,,fe,,,,,,-) = -7(67 - /,) + + 1 In [/(O, yM)] , (77) 

- n^oo pn - p p 

and the unique global maximum of |/| at (0, 2/^(0)) on the curve z{s) excludes any phase transition for this 
antiferromagnetic model. 

Finally, notice that Lemma |H1 shows that the only extreme point of |/p on the y-axis is at the y- 
intercept of the curve z{s). This extreme point must be a saddle point of the holomorphic function /(-z). 
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hence {0,yz{0)) is a global minimum of v{y) = f{0,y). An explicit formula for the Cournot free energy 
can then be written 

1 ln(q — q) 

F{P, a, b, c, g, g) = - - (g + q) [b{q + g) - 2(a - c)] + ' 

1 . r .^ nx (78) 



+ — mm 

P ?/£(— co,oo) 



„ / r _ (q-q) _ ]\ 

(3by + In ( sinhc Pb{q — q)y + P — ^^'^ {b{q + q) — (a — c)} j 



where y was translated and rescaled (this does not alter the result since y ranges over all real numbers) to 
get rid of square roots, and sinhc(M) = smh.{u)/u with sinhc(O) = 1. 

Now for a final note on the expected value of the variable g^. The minimum value in (j78|) occurs when the 
derivative of the smooth function within the minimum is zero, which is at ym = 1/2(0). Since the minimum 
is a global minimum, the second derivative with respect to y is nonzero. By the implicit function theorem, 
the minimum point ym is locally a smooth function oi h = a — c. As a result, the partition function in the 
form (f77j) is seen to be a smooth function of h. 

The function Fn{h) in ()53p is convex in h via Holder's inequality, and as a result 

F„ih)-F„ih-S) ^ ^, ^ F^ih + S)-FM ^73) 

6 — n — ^ 

where 5 > Q. Since F{h) is a smooth function of h, taking the limit of (|7^ as n ^ 00 results in 
F'{h) = lim„_ooF^(/i). Using §^ 

F'M=l-v/b-i{i)''J{2pVb), (80) 

where (■)^ is a (one-dimensional) Fourier transform of the Gibbs measure for n agents, (/)„ = (/ / e'^^)/2„, 
which is generated by (|52j) . The measure (■)^ is generated by the form of the partition function in (j59p . 
Taking the limit of (|8(Jj) as n — > 00, 

F'{h) = h/{2b) + lim {2t/{QVb))^/{2l3y/b). (81) 

n^oo 

The latter limit can be evaluated by deforming the integral in the complex plane as before (change t 
to z), and using Lemma IHl The sequence (■)^ of measures converges to a 'delta function' or evaluation 
measure at (0, ?/„,). Such measures are precisely the characters, meaning they are homomorphisms on the 
'algebra of observables' and as such are extreme points. This shows that there is no phase transition in 
this antiferromagnetic model. It is worthwhile to elaborate on such technicalities to show that the weak 
convergence of the Fourier transformed measures implies that the sequence of measures (■)„ converges 
weakly. With LemmaEl this implies the Gibbs measures (■)„ converge weakly to the infinite-agent measure 

(■)• 

Lemma 9. The weak convergence of the sequence (•)^ of measures implies the weak convergence of the 
sequence of Gibbs measures (■)„. 



Proof. Let 2t„ be the algebra of observables for n agents, which is simply L°°{\q, q\^, dqi ■ dqn). The algebra 
of observables (functions) for an infinite number of agents, 2too, is the inductive limit (von Neumann algebra 
sense) of the finite-agent algebras 2t„. As a result, the 2l„ can be considered subalgebras of 2loo and any 
weak limit of finite-agent measures is defined on 2too. 

Consider two distinct limit points (■), (■)' G 2t^ which are limits of subsequences (■)„i and (■)n2! respec- 
tively. There is a function g G 2loo on which (■) and (■)' differ. By approximating g closely enough with a 
function in some 2t„Q we may assume g G 

""^■^This is the one-dimcnsional Fourier transform of the Gibbs measures (•)« in the variable s„ = (gi + • • • + (7ri)/v^ that 
is generated by the partition function in l)59|l 



25 



The integral {g)n, over the g^, 1 < j < rii can be manipulated into an integral over the variable i as was 
done with ()59|). upon which it becomes an integral {G)^_, where G = G{t) is a function independent of n^. 
This shows {g) = lim„^(G'),^^ and {g)' = limn2{G)^^, which implies the assertion. □ 

Remark. With q = (gi, . . . g = g{q), note that (j^UI) inserts the Fourier transform ^ (of the Gibbs 
measure) into the {g)n integral, and ()57|) shows what the dual of the Fourier transform, ^* does to g. 
The Fourier transform of the Gibbs measures need only converge on the subalgebra of functions generated 
by the ^* [21^] in order for the sequence of Gibbs measures to converge. 



The expected value {qi) is then 



From (|7H|l . we can see that 

which simply states that for completely irrational behavior (i.e., 'high temperature'), the agents will act 
randomly, and the Gibbs-Cournot measure will be uniform. As such, the expected value (g,) over the 
interval [q, q] is the average 7 = {q + q)/2. 



In the completely rational limit (i.e., 'zero temperature'), if Q/2+7 — /i/(26) > and ?7 < (equivalently, 
if 7 < h/{2h) < q), then the curve g{y) which equals the right side of (f7B|) has a positive limit as ?/ — > —00. 
In this case, the line h{y) = y/f3 intersects g{y) at a point with y-value ym which has a finite limit as 
P ^ 00. Thus ym/ /5 — > as /5 — ^ 00 and 

lim F'{h) = 4 for 7 < h/{2b) < q. (84) 

In the case when h/{2h) > q, the curve g{y) is always negative-valued. As a result, ym —00 as /5 ^ 00. 
In this case, g{y) can be used to evaluate ym/ P and 



l„„nft) = 4+ (1 + 0^-1;)=? for V(2(,) > 



Since the integral in is even in r/ = 67 — h/2, lim^^oo -^'(^) = Q when h/{2b) < q. If 67 — h/2 > 0, 
then g{y) in ([7^1) is shifted up and g{0) > 0. With this, g{y) can only intersect the line y/P at positive ym- 
It is then seen that ym > when h/ (2b) < 7 and ?/m ^ when h/{2b) > 7. 



In summary, we see that the expected value (qi) in the irrational case is the average value {q+q)/2. In the 
completely rational case, we recover the classical Nash equilibrium value (qi) = h/ (26) when q < h/ (26) < q. 
When h/ (26) < q {> q) then (g^) will be g (g). For finite /3, the equilibrium (g^) will be smaller (larger) than 
h/(2b) if h/{2b) itself is larger (smaller) than the average 7 = (g + g)/2. Furthermore the graphical analysis 
below (fTHj) shows that \ym\/P is a decreasing function of /5, and ym < 0(> 0) when h/{2b) > 7(< 7). As 
a result, the agents deviate farther from the Nash equilibrium h/{2b) and closer to the average 7 as /? 
decreases (i.e., as deviations from rationality increase). 



"'^''in particular, ^*[g]{Vbt) = J cxp(i\/6ts„)£^[(7(s„, s;|;)Xq(s, s;|;)|s„]c?s„, where s:^ are integration variables orthogonal 
to s„, xq is the indicator function on q E [— Q/2, (5/2]", and i?[-|s„] is the conditional expectation over the subalgebra of 
functions in the variables s;!-. 
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The volatility (i.e., susceptibility) 

1 " 

X = -F"{h) = lim - (86) 



/3 n^oo n 



can be evaluated in a similar manner. Let 



G{K 1/) = ^ (- coth 1/ + ^) + Q (^7 - I) - I (87) 



which equals zero at points {h,ym{h)). Then 

1 dym Q 



(3 dh /?6Q2(sinh-2^^_^^2)_2' ' 

In the irrational limit, lim^^o2/m = and lim^^o = (sinh^^y^ — U^) = —1/3. Differentiating (jH^ results 
in 

Q^(sinh-^y^-y^2)/2 

and as a result 

limx=— . (90) 

In the completely rational limit, either j/^ converges to a finite value (when (g + g)/2 < h/{2b) < q) or 
goes to — oo (when h/{2b) > q). If converges to a finite limit, then the denominator of (j89p goes to 
-oo and lim/3^00 X = 0. If ym goes to -cx), then x < ((5^/4)| sinh"^ - y-^\ and lim/3_ooX = 0. Using 
that ()59|) in even in 77 as before, 

lim X = 0. (91) 

In summary, when agents behave randomly in the irrational limit, the volatility is simply the variance 
of a uniform random variable. Each agent becomes decorrolated from the other agents and the volatility 
is the limiting-average of the volatility of each agent. When agents are completely rational they all output 
at the Nash equilibrium values and do not deviate from that, which is reflected by zero volatility. 
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